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Abstract 

1^ . In this paper we calculate the 0{as) corrections to the xp- and y-distributions of 

lepton pairs produced in collisions of longitudinally polarized hadrons. The numer- 
ical importance of these corrections is studied and consequences for the extraction 
^ ■ of the polarized sea quark distributions from a measurement of the longitudinally 

H . polarized Drell- Yan cross section are discussed. 
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1 Introduction 

The production of lepton pairs in hadron collisions, the Drell-Yan process 0, is one of the 
most powerful tools to probe the structure of hadrons. Its parton model interpretation is 
straightforward - the process is induced by the annihilation of a quark-antiquark pair into 
a virtual photon which subsequently decays into a lepton pair. The Drell-Yan process 
in proton-proton or proton-nucleus collisions therefore provides a direct probe of the 
antiquark densities in protons and nuclei. 

The naive quark parton model predicts scaling of the Drell-Yan cross section in 
the variable r = M'^/S, where M is the invariant mass of the lepton pair and a/S* 
is the hadron-hadron centre-of-mass energy. Apart from the invariant mass distribu- 
tion dcr/dM^ one usually studies the distribution of the lepton pairs as function of the 
Feynman-parameter 

-. = ^ (1) 

or of the hadron-hadron centre-of-mass rapidity 

2 qo-Qz 

where q denotes the four-momentum of the lepton pair in the hadron-hadron centre-of- 
mass system. The resulting distributions 

dcr dcr 

and 



dMMxp dM^dy 

can be directly related to the x-dependence of the parton distribution functions in beam 
and target. Moreover, most fixed target experiments have only a limited kinematic cover- 
age in xp or y, such that only these distributions can be measured without extrapolation 
into experimentally inaccessible regions. 

A summary of the experimental data on the Drell-Yan process in unpolarized pp-, 
pN-, pN- and vrA^-collisions can be found in 0. These data are used to determine the 
large-x behaviour of antiquark distributions in the proton (e.g. 0, Q) and to study the 
parton structure of the pion 0. In recent times, the difference between Drell-Yan cross 
sections in pp and pd collisions [|] has provided valuable information on the asymmetry 
between the u- and J-quark distributions in the proton J^. 

Up to now, the only experimental information on the polarized parton distributions in 
the proton comes from deep inelastic lepton-nucleon scattering experiments |p, measuring 
the polarized structure function (yfi(x, Q^). This structure function probes one particular 
combination of the polarized quark distributions, its mere knowledge is therefore insuffi- 
cient for a determination of all different quark and antiquark distributions. This lack of 
knowledge is reflected in the results of recent fits of polarized parton distribution func- 
tions p, 0, |Tl[| . While the polarized valence quark distributions can be determined with 



some precision from the experimental data, only little information can be extracted on 
the total magnitude of the polarized sea quark distribution; the fiavour structure of the 
polarized quark sea is at present completely undetermined. 



Experimental information on the polarized antiquark distributions can be obtained 
from several observables. Charged current exchange in polarized electron-proton scatter- 
ing can probe the polarized weak structure functions |T^, which are sensitive to a different 



combination of quark distributions than their electromagnetic counterparts. Combining 
electromagnetic and weak structure functions in a global parton distribution fit, it should 
in principle be possible to disentangle the different quark and antiquark flavours. A direct 
probe of the polarized sea quark distributions would be the Drell-Yan process in polarized 
proton-nucleon collisions |T^. Such a measurement could be feasible at RHIC [|T^] or at 



HERA p|. 



The rather large QCD corrections to the unpolarized Drell-Yan process suggest, that 
a reliable interpretation of the Drell-Yan cross section in terms of partonic distribution 
functions is only possible if QCD corrections are included. This paper aims to derive the 
0{as) corrections to the xp- and ^-distributions of lepton pairs in longitudinally polarized 
hadron-hadron collisions and to study the impact of these corrections on the Drell-Yan 
cross section. 

For the unpolarized Drell-Yan process, QCD corrections have been calculated up to 
0{a1) for the invariant mass distribution |TB|, ^ and up to 0{as) for the rapidity and 



XiT-distributions [|T8|, [T^. For the latter, parts of the 0{a1) corrections are known |2i 
Inclusion of the (9(Q;s)-corrections changes the Drell-Yan cross section at fixed target 
energies by up to 50%, the C(a^) corrections amount up to 10%. At small transverse 
momenta of the lepton pair, the QCD corrections are dominated by the emission of soft 
gluons. The effects of soft gluon emission at arbitrary order in perturbation theory have 
been studied for the Drell-Yan process in [^ , yielding a resummation of the dominant 



corrections to all orders. 

All these results can in principle be applied to the production of vector bosons at 
hadron colliders. However, in the case of W^-boson production, one has to take into 
account that PF-bosons in hadronic collisions can be identified only in their Ivi decay 
mode, where just the lepton is observed. All above calculations cannot make predictions 
for the kinematic distributions of the vector boson decay products in the hadron-hadron 
centre-of-mass frame. Separate calculations of the next-to-leading order QCD corrections 
to the transverse mass distribution of the W^-boson [p2|, to the polarization density matrix 



of its decay products |Q and to the production rates of 1^-1-0,1 jets |^ in unpolarized 
hadronic collisions have been performed. 

QCD corrections to the Drell-Yan process in polarized hadron-hadron collisions have 



been derived up to 0{as) for the invariant mass distribution in longitudinally [53, EB 



and transversely [^ polarized collisions. The 0{as) corrections to the rapidity distribu- 



tions are known for longitudinally [^ and transversely |28[ polarized collisions. Further- 
more, a resummation of the effects of soft gluons in longitudinally polarized collisions has 
been performed in |26|| . 

A fully consistent study of the Drell-Yan process at next-to-leading order was until 
now only possible in the unpolarized case, as the polarized parton distributions could 
only be determined at leading accuracy. With the recently calculated polarized two-loop 



splitting functions [^, the polarized distributions can now be determined to next-to- 



leading order from fits [^ [TD|, |n[] to polarized structure function data. Having a complete 



calculation of the two-loop splitting functions available, it is now furthermore possible 



to define consistent scheme transformation prescriptions p9| , |TT| for parton distributions, 
splitting functions and parton level cross sections at next-to-leading order. The scheme 
dependence of the parton level cross sections for polarized deep inelastic scattering has 



lead to some discussion [^] about the partonic interpretation of measurements @] of 
polarized lepton-nucleon scattering. The scheme dependence of the QCD corrections to 
the polarized Drell-Yan process has been addressed in ||2^, ^ . 

This paper is organized as follows: in Section |^, we derive the (9(as)-corrections to the 
XiT-distribution of Drell-Yan lepton pairs produced in longitudinally polarized hadron- 
hadron collisions. We discuss the representation of 75 used in our calculation and motivate 
our choice of factorization scheme. Furthermore, a brief rederivation of the corrections to 
the invariant mass and ^-distributions is given. We numerically estimate the importance of 
these corrections at fixed target energies in Section ^. Section ^ contains a brief summary 
of our results and concluding remarks. In the Appendix, we collect several identities used 
in the calculation. 



2 Perturbative corrections to the polarized Drell-Yan 
process 

The QCD corrections to the XiT'-distributions in the unpolarized Drell-Yan process were 
originally derived using dimensional regularization by Altarelli, Ellis and Martinelli |l^ . 



In a later publication, Kubar et al. |T9| confirmed the results of |T^ and derived the QCD 
corrections to the y- distributions. The calculation of Kubar et al. was performed using a 
non-zero gluon mass as regulator. The results obtained for the xjT'-distributions with this 



method were considerably simpler than the original results of [Iq]. In the calculation of 



the QCD corrections to the xp- and ^/-distributions in the polarized Drell-Yan process. 



we shall work in dimensional regularization and follow closely the method of [|I^. Using 
several identities listed in the Appendix, we are able to simplify our results into a form 
similar to the unpolarized results of [|19[. In the following, we give an explicit derivation 
of the corrections to the Xi?-distribution. Furthermore, we briefly summarize the results 
of a rederivation of QCD corrections to the y-distribution within the same method. These 



corrections were originally obtained by Weber p6[ as a 'by-product' in the calculation of 
the soft gluon resummation in the polarized Drell-Yan process. 

Disregarding correlations between the final state lepton plane and the direction of 
the incoming hadrons, it is sufficient to consider the production of an off-shell photon of 
invariant mass M^, the polarization sum of the outgoing photon can be taken to be —g'^'^. 
The lepton pair production cross section can be retrieved by multiplying the results with 
a factor a/{3TcM^)dM'^ for the decay 7* -^ l+l' . 

Only the lowest order contribution to the Drell-Yan process (Fig. |l|(a)) yields a finite 
cross section. All higher order parton level cross sections contain infinities. To make these 
infinities explicit, we work in dimensional regularization with d = 4 — 2e dimensions. This 
treatment faces the problem of describing 75 in n 7^ 4. We shall use the 75-prescription [^ 
of 't Hooft, Veltman, Breitenlohner and Maison (HVBM), which has been consistently 



used in the derivation of the next-to-leading order corrections to the polarized splitting 
functions ||2^. In this prescription, one restricts the anticomniutation property of 75 



to the physical four dimensions, while 75 commutes in the remaining n — 4 dimensions. 
The major drawback of this formalism in the MS-scheme is the non-conservation of the 
flavour non-singlet axial vector currents [^ due to a non-vanishing first moment of the 
corresponding non-singlet NLO splitting function AP^g _|_. To restore the conservation of 
this non-singlet axial vector current, a further scheme transformation ^^ of the results 



obtained in the HVBM formalism is needed P^jp], We have to apply the same scheme 
transformation to our results for the parton-level Drell-Yan cross sections at 0{as)- 
We define the hadron cross sections: 

dAa = ]- (d(7++ - d(j+-) , da = ^ (da++ + d(7+-) , (3) 

where (+) and (— ) denote positive and negative hadron helicities. In the parton model 
picture, these cross sections can be expressed as convolutions of (polarized) parton distri- 
butions [A]/j with parton subprocess cross sections [A]o-jj: 

d[A]cr = / dxidx2X]d[A]aij(xi,X2,/u|) [A]/i(xi,/i|,) [A]/,(x2,/i|.). (4) 

We have defined 

dA(T = - (da++ - da+-) , da = - (da++ + da+-) , (5) 

where (+) and (— ) now denote positive and negative parton helicities. 

The lowest order bare cross section for the qq -^ 7* contribution (Fig. 0.(a)) reads in 
d = 4 — 2e dimensions: 

Aa° = -^^(1 + 6)5(1 -z), (6) 

where we define z = M'^/s, with y/s being the parton-parton centre-of-mass energy. The 
factor (1 + e) is common to all contributions and forms part of the overall normalization 
of the cross section. 

The virtual gluon correction (Fig. |lKb)) to this process yields 

^.y 47r2ae2 « f4:nij,^Y 1 / 2 3 „ gA 



3M2 ' ' ' ' 271 ' \ M^ J T{1 - e) \ e^ e 

The correction due to real gluon emission (Fig. |l].(c,d)) takes the form: 



3M2 ' ' 27r V M2 / Fd 



e 



X 



/' dz^zni - z,)-^ [(1 + 3e) f ^^ + -^^) 
Jo V ^1 1 — 2i / 



2z 



Zx l-ZiJ [I- zYzi{l-Zi, 



^This scheme transformation was first considered by Weber [ p6| to restore the conservation of quark 
hehcity in the gg-contribution to the polarized Drell-Yan process. 





^y^ 





(e) 




Figure 1: Parton level subprocesses contributing to the Drell-Yan process up to 0{as) 
lowest order gg-anniliilation (a); virtual (b) and real (c,d) gluon corrections to qq- 
annihilation; quark-gluon Compton scattering (e,f). 



where we have replaced the parton-parton scattering angle O* by zi = {1 + cosO*)/2. 
Terms proportional to e have been omitted if they do not yield a finite contribution to the 
final result. Finally, the contribution from the quark-gluon Compton scattering process 
(Fig. |l|.(e,f)) reads: 



2 ^,^2 



Aa^'^ 



Arc^ae 
3M2 



l(l + e)^T^''"^ 



27r 



M2 / Fd-e^ 



z^+Hl-z' 



,l-2e 



X / dziz^ '(1 - zi] 
Jo 



-^— + 2 - (1 + e)- - (2 - zi)(l - ^) 
,l~z)zi Zi 



(9) 



Using eqs. ( |A.1| ), one can make the poles in e in the annihilation and Compton contribu- 
tions explicit: 



Aa 



M = Aa^'^ + Aa^'^ 

3M2 ^ ^ ' 2t:^\ M2 ) F(l-e) 
X f^ dzi J i^-^5{l- z) (5(zi) + 5(1 - zi)) 



+5{l-z) 



1/1 1 

+ 



e V(^i)+ (l-^i)- 



^/ln^\ ^^ln(l-;.i^ 



Zl 



l-zi 



Inzi ln(l — zi] 



1-^1 



2:1 



+ [5iz^) + Sil-z^)] 



+ 



1 



In- 






+ 



(zi)+ il-zi)+ 



.2^,„2 



A(T 



1,C 



47r ae- a, 

" 3M2 ^ ^^ 27r ^ I M2 



-1-2 + 



z 
2 



1 



\-z 
2(1-^)1 



Xtvz 



(10) 



rfi 



X / Azxzl - -5(2;i) [22; - 1] + (5(zi) 







(22-1) In. 



2-22 



l^i) + 



[22-l]+2(l-2)-(2-2i)(l-2)^ 



;ii) 



The XF-differential cross section can now be obtained by multiplying the above expressions 

by 

/ &Xp S (xp - [zi (xi + Z2) (1 - 2) + ZXi - X2]) . 
Using the definitions 



.0 _ ™0 

"1 -^2 



Xp — xX — X9 and t — x^x 



0™0 
1"''2' 



one can fix 



^^1 = ■^{^F + \Jx\ + 4r , , 



X 



2 - 2 I ~^P 



x'i + Ar 



12) 



for the XiT'-distributions. The parton level cross sections Ao"^'"^ and Ao"^'*" can then be 
rewritten using eqs. ( |A.3| )-( [70D . 

The resulting parton level cross sections still contain divergences associated with 
collinear singularities in the initial state. These collinear singularities are process indepen- 
dent and are absorbed into the bare parton distribution functions by a mass factorization 
procedure. The only distribution relevant in our case is the polarized quark distribution. 
Performing the mass factorization in the MS-scheme and making a scheme transforma- 
tion p9| , |2^ to remove spurious terms due to the HVBM prescription for 75 yields for the 
bare quark distribution up to 0{as): 



a, /■! dx' 



Aq{x) = Aq{x,fil)--^ -j- 

ZTT Jx X 



1 1 



with 33 



AP,,(x) 



AP,,{x) 



as /■! dx' 
'2^Jx ~ 



C 



1 1 



er(l-e 

,2\ <= 



477/^2 



^1 AP,,(a:') 

PF / 



AP,,(S') + A2gg(x') 

AG{x/x'), 



Aq{x/x') 
(13) 



(1-x). 
Tp [2x - 1] . 



1 - X + -S{1 - x) 



AZgg{x) 



-4Cp(l - x) 



After mass factorization, the polarized cross sections for the Drell-Yan production of 
lepton pairs at NLO can be expressed in a form similar to the unpolarized cross sections P] : 



dAa 
dM^dxF 



X- 



9APS 



E' 



da;i / da;2 



dAa, 



(0) 



a;i,X2) + 



as dAa^^^ 



dA'PdxF 27rdM^dxF 

+Aqi{xi,ijl)Aqi{x2, Ai|)} 



+ 



as dAgW 



M2 






\^Aqi{xi, idp)Aqi{x2, IJ.I) 



27r dM^dx I ^i' ^2, ^- I AG(xi, /i|) {Agi(a;2, /i|) + Agi(x2, nl)} 



+ (1 



(14) 



with 



dAa, 



(0) 



dMMx^ 



(1) 
gg 



(a;i,a;2) 



M2 



dAa^^ 






(5(a;i — xl) 6{x2 — x^) 

J/1 n^ Jjn 



da. 



(0) 

gg 



dMMxi 



-(a;i,a;2) 



(15) 



-C7, 
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1 "*" 2 
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8 + 2Li2(a;°) + 2Li2( 



x°) 



+ 



+ ln2(l 
'5{xi - 
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+ ln2(l -x°) + 21n 
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1 
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1 "I" 2 



2 
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(16) 



(1) 
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dM2dXir 
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Xly X2i T~ 



-T, 



(Xi + X2)Xi 



+2(xi 
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:i7) 



where 



AG'^(xi,X2) 

AH^{xi,X2) 



{xi + X2) {xi ^^2 ^ + x^a;2 

xfxKx? + X2)(Xi + X2) 

2 



a;ia;2(a;i +0:2) 
2a:?X2 — X1X2 



XiX2(xi +0:2) 

a;i(x? + 0:2) (2;2 - ^2) + 2x?X2(a;i + X2) 



Inserting 



Jdys{y-l 



In h In 

X2 



xfxKxi + X2y 

Z + Zi — ZZ\ 



into eqs. (|lO|)- ([Tl|) and using identities (|A.5| ),( |A7fD - 

corrections to the rapidity distribution of the Drell-Yan pairs. The variables 



x\ 



v^e^ 



X2 



1 — Zl + ZZ\ . 

one can derive the 0{a^- 



(18) 



VT-e" 



are in this case defined by 



2 x% 



and 



0^0 



X 1 Jj 



\-^2- 



They are in general different from the x\ and x\ occuring in the Xir-distributions. 



These corrections have already been calculated in |2g]; we find numerical agreement 
with the results (keeping in mind that a different factorization prescription for the gluonic 
contribution was used in p6|), our expressions are however considerably simpler: 
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(19) 



with 



"^^^ ix X) 



dAa, 



(1) 

99 



where 
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As mentioned above, the (9(Q;s)-corrections to the xp- and ^-distributions in the unpo- 



larized Drell-Yan process were derived in [jT8|, 0. A summary of these formulae in the 
MS-scheme in our notation can be found in [^|. 

Finally, integration of eqs.(|TO|)-(|lT]) with respect to Zi reproduces after mass factor- 



ization using eq. ([T3| ) the correction terms to the Drell-Yan mass distribution [g^, ^ |27 
These read in the MS-scheme: 
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3 Numerical results 

In the following, we will illustrate the impact of the QCD corrections to the polarized 
Drell-Yan process derived above by numerically evaluating (14) and (19) for recent pa- 
rameterizations of the polarized parton distribution functions. 

When using the polarized GS parton distribution functions 0, we take A„ 4 = 
231 MeV, the corresponding unpolarized cross sections are then evaluated using the un- 
polarized MRS parton distribution functions set A' 0]. The polarized GRSV distribu- 
tions [jlOl are consistently used in combination with the unpolarized distributions from 
GRV 0] and for A„ 4 = 200 MeV. If not stated otherwise, we shall always use fip = M; 
the strong coupling constant a^ is evaluated at fip- 



^ As the normalization of the gluon polarization sum used in |5|, y_8|; |l9|] deviates from the normalization 
conventionally used in dimensional regularization, eqs.(A8) and (A20) of [^ have to be modified according 
to |]34| to conform to the normalization used in the NLO splitting functions and parton distributions. 
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Figure 2: The Drell-Yan asymmetry A{M) in pp collisions with Ebeam = 
different parameterizations of the polarized parton distribution functions. 



820 GeV for 



All results in this section are obtained for ^/s = 39.22 GeV, corresponding to a fixed 
target experiment in the HERA proton beam [^ with Ef,eam = 820 GeV. 



At present, all polarized parton distributions are fitted to data on the polarized struc- 
ture function gi{x,Q'^) only. As this structure function is dominated by valence quark 
contributions for x > 0.1, the knowledge on the polarized sea quark distributions in the 
large-x region is very poor. This uncertainty is illustrated in Fig. |^, which shows the 
next-to- leading order predictions (Eq. (ESI)) for the Drell-Yan asymmetry 



A{M) 



dA 



a 



da 



dM / dM 



obtained using the polarized parton distribution functions of [^ (GS(A-C)) and 
(GRSVs,v). At present, neither magnitude nor sign of this asymmetry can be predicted. 
It should therefore be clear, that all numerical results for the polarized Drell-Yan process 
presented in the remainder of this section do not make solid predictions for observable 
asymmetries. The purpose of these numerical calculations is only to illustrate the magni- 
tude of the next-to-leading order corrections and to demonstrate the perturbative stability 
of the predictions. 

In the following, we shall always study the cross section asymmetries 



A{xf) 



dA(T 



da 



dMdxp / dMdxF 



and 



My) 



dAa 



da 



dMdy / dMdy 



As the unpolarized parton distributions are known to some level of accuracy, it is pos- 
sible to predict the cross sections in the denominator of A[xf) and A{y). These cross 
sections are shown for proton-proton collisions {y/s = 39.22 GeV) in Fig. ^, which is ob- 
tained with the MRS (A') parton distributions [Q. The prediction using the GRV parton 
distributions [Q agrees within 3% over the whole kinematic range. 
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Figure 3: The unpolarized Drell-Yan cross section in pp collisions. Curves are for different 
invariant masses M; solid: 4 GeV, dashed: 6 GeV, dotted: 8 GeV, dot-dashed: 12 GeV. 
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Figure 4: The Drell-Yan asymmetry in pp collisions for the polarized GS(A) parton 
distributions. Curves are for different invariant masses M; solid: 4 GeV, dashed: 6 GeV, 
dotted: 8 GeV, dot-dashed: 12 GeV. 



It can be seen from Fig. |] that the total Drell-Yan asymmetry A{M) changes with 
the invariant mass M. The change is however not uniform over the whole range in 
xp or y as illustrated in Fig. ^ for the polarized GS(A) parton distributions. The xp- 
dependence of the asymmetry can be understood from eq. ([l^). The asymmetry around 
xp = always probes the polarized parton distributions in beam (1) and target (2) around 
T. If xp is increased towards its positive kinematic limit, x^ increases towards 
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Figure 5: Contributions of the individual parton level subprocesses to the polarized Drell- 
Yan cross section. The subprocess cross sections (short-dashed, dot-dashed and dotted 
lines) and their sum (solid line) are normalized to the full unpolarized cross section at 
NLO. The long-dashed line shows the asymmetry evaluated at leading order. 
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Figure 6: Sensitivity of the Drell-Yan asymmetry on variations of the factorization scale. 



1 while x\ decreases only slightly. Decreasing xp towards the negative kinematic limit 
yields a small decrease of x\ while x\ increases towards 1. The region xp > ^ therefore 
probes the x-dependence of the parton distributions in the beam hadron, while Xi? < is 
sensitive on the target hadron structure. Moreover, different invariant masses correspond 
to different x intervals probed. A similar argumentation, based on eq. (|18|), applies to the 
y-dependence of the asymmetry. 

The contributions of the individual subprocesses (gg-annihilation at leading and next- 
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to-leading order and quark-gluon Compton scattering) to the polarized Drell-Yan cross 
section are illustrated in Fig. ^. All curves are obtained with the polarized GS(A) parton 
densities and are shown for M = 8 GeV. All polarized subprocess contributions are 
normalized to the full unpolarized cross section at next-to-leading order. The relative 
magnitude of the individual contributions is similar to the unpolarized Drell-Yan process. 
The 0{as) correction to the gg-annihilation process enhances significantly the lowest order 
prediction while the quark-gluon Compton process contributes with a sign opposite to 
the annihilation process. However, the relative magnitude of annihilation and Compton 
process depends on the magnitude of the gluon distribution at large x, which is completely 
undetermined at present [^ |lO| . This uncertainty prevents a sensible prediction of a K~ 
factor between the cross sections at leading and next-to-leading order. 

Furthermore, Fig. |^ shows the prediction for the Drell-Yan asymmetry if both polar- 
ized and unpolarized cross section are evaluated using the leading order qq, (O)-coefficient 
function only (long-dashed line). The sizable shift between the leading and next-to- 
leading order predictions for the Drell-Yan asymmetry seems to contradict the claim [^ 
that the next-to-leading order corrections, although large in polarized and unpolarized 
cross sections, would largely cancel in the Drell-Yan asymmetry. 

The change of physical quantities such as cross sections or asymmetries under varia- 
tions of the (unphysical) mass factorization scale fip enables an estimate of the numerical 
importance of unknown higher order corrections. We quantify this uncertainty in Fig. |^, 
showing the Drell-Yan asymmetry for hf = 0.5, 1, 2M. It can be seen that the abso- 
lute value of the asymmetry changes less than 0.01 in the central region. Even in the 
region where xp (or y) approaches its kinematic limit, the absolute value of the asym- 
metry changes only by a maximum of 0.015. This stability is crucial for an estimate of 
the theoretical error on the polarized sea quark distributions extracted from the polarized 
Drell-Yan process. 

A final point briefly concerns phenomenological aspects of the polarized Drell-Yan 
process, although this is not the principal aim of the present paper. Figures |^ and ^ show 
the Drell-Yan asymmetry in polarized proton-proton and proton-deuterium collisions for 
M = 8 GeV as predicted by using the different polarized parton distributions suggested 
in [^ and ||T0|. It can clearly be seen that the different predictions vary within a band 
of ±0.06 in pp and ±0.04 in pd collisions. This uncertainty reflects the present lack 
of knowledge on the polarization of the light quark sea in the proton at x ^ 0.1. The 
difference between the predictions is furthermore substantially larger than the theoretical 
uncertainty due to the choice of the factorization scale. 

Due to the smaller overall valence quark polarization in deuterium, one finds the 
asymmetry in pd collisions to be systematically smaller than in pp collisions. Although 
the normalization of the asymmetry in the proton direction [xp, y > 0) and in the deuteron 
direction [xp, y < 0) is slightly different, its shape is still very similar in both regions. 
This reflects the assumption of an SU (2) /-symmetry of the polarized antiquark sea, which 
is used in all polarized parton distributions displayed in the figure. A measurement of 
the Xf-{ot I/-) dependence of the Drell-Yan asymmetry in pd collisions could therefore test 
this assumption. 
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Figure 7: Predictions for the Drell-Yan asymmetry at M = 8 GeV in polarized pp 
collisions {Eheam = 820 GeV) obtained for different parameterizations of the polarized 
parton densities. 
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Figure 8: Predictions for the Drell-Yan asymmetry at M = 8 GeV in polarized "pd 
collisions [Eheam = 820 GeV) obtained for different parameterizations of the polarized 
parton densities. The proton is taken in xi-direction. 



4 Conclusions 

We have presented a complete calculation of the 0{as) corrections to the xp- and y- 
dependence of the longitudinally polarized Drell-Yan cross section in the MS-scheme. 
The results for the ?/-dependence and the mass distribution agree (after suitable change 
of the factorization scheme) with earlier results in the literature [^, |2^, ^. We have 
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demonstrated that the impact of these corrections on the polarized Drell-Yan cross section 
is very similar to the impact of the unpolarized corrections. The corrections due to the 
individual subprocesses turn out to be sizable. However, a K-iactoi between leading 
and next-to-leading order results can at present not be sensibly predicted due to the 
uncertainty on the magnitude of the polarized gluon distribution. 

The corrections derived in this paper are of particular importance for the extraction 
of the polarized sea quark distributions from measurements of the Drell-Yan process. 
As discussed above, a measurement of the asymmetry as function of xp or y can be 
directly inferred into the x-dependence of the polarized sea quark distributions. The 
results obtained in this paper can be furthermore immediately - keeping in mind the 
limitations discussed in the introduction - applied to the production of massive vector 
bosons in polarized proton-proton collisions at RHIC . 



With the knowledge of the next-to-leading order corrections, it is furthermore possible 
to quantify the uncertainty of the theoretical prediction and hence the theoretical error on 
a measurement of the polarized sea quark distributions. We have demonstrated that the 
absolute value of the Drell-Yan asymmetry in the central region varies by less than 0.01 
under variation of the mass factorization scale, demonstrating the perturbative stability 
of the prediction. This uncertainty is far smaller than the deviation between predictions 
obtained for different parameterizations of the polarized sea quark distributions. 

A Identities for "+"— functions 

The e-expansion of the singular terms in the parton level cross sections of Section |^ yields 

Integration over zi yields the following "-I-" -functions, all defined on the interval [0; 1]: 

1 1 




) 



1 



and (5-functions: 

Furthermore, various products of the above terms appear. By construction, we have 
z = {x1x2)/{xiX2) while zl is fixed by the definition of x^? or y respectively. As the parton 
level cross section is always integrated over the parton momentum fractions Xi and X2, it 
is desirable to re-express the above functions in terms of 

6{xi-x[), 6{X2-X2) 
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and 



ln(l — x^/xi) 



with 
with 



dxi- 



f{Xi 






'a^.lM^lMl, 



tXj 1 tXi A 



tXj 7 tXy A 



•d.,/wf'-'<^-^°/^->' 



X 7 X „■ 



= / da;,- 



f{xd - f{x1 



ln(l — x^Xi) 



/y . lyyi 



(A.2) 



[{Xi -X'i){x2-X^)] + 



with 



da^i / dx 



f{Xl,X2) 



4 [(a;i-x?)(x2-a;^)] + 



^ dx C dx ^'^^^'^^^ ^ •^'^^1' ^2) - /(a^i> ^2) + i\A, xl) 



{xi -x\){x2 -xl) 
We first consider the Xi?-distributions. All terms containing 5{1 ~ z) can be simplified 



by 



(5(1 — z) 5 {xp — [zi{xi + 3:2) (1 ~ z) + ZXi — X2]) 






„ -fr(5(a;i -a;i)5(a;2 -^2) (A. 3) 

1 "I" 2 



and subsequent integration over zi. For the remaining terms, we integrate over Zi using 
the XiT-fixing (5-function, which yields 



x" + X2 



(1 - 2;)x2(a;i + ^2) 



{X2 - X2} 



and 



1 



iAj "I I iX/Q 



(1 - z)xi{xi + 3:2) 



I iX/ ]^ iX/ -| 



As the integrations over the parton momenta xi and X2 are independent from each 
other, we can perform them in arbitrary order. For terms containing S{zl) and l/{zl)^ 
only, we choose to integrate X2 first. For terms containing 6{1 — z*) and 1/(1 — z*)-|- only, 
Xi is integrated first. We can then make the following replacements: 



6izl) 
1 

6{1 - zl) 

1 

(1 - ^1*)+ 



(1 ~ z)xl{xi +x^) 



(1- 


-Z)X2{XI+X2) 
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-Z)X^,{X',+X2) 


(1- 


X^, + X^2 
- z)xi{Xi+X2) 



6{X2 -^2), 



[X2 — 3^2) + 
5{xi -X?), 



+ 5(a;2 — ^2) In 



0\ i„ l-^i "T •^2)1-'- X2JX1 



a:2(a;i + a;2)(xi — xi) 



(A.4) 



Xi + X2 



1 I r/ _ On 1 (^1 + ^2)(1 ~ X]^)X2 

(Xi-x5)+ "^ X?(x5 +X2)(X2 -X^) 



The above identities for 5{zl) and 5(1 — z^) can be safely applied to terms containing 
"+"-functions in [1 — z). These can then be rewritten as 



(5(xi — x\) 



5{xi — x\) 



X2 



tLn 



[X2 — X2) + 



+ X2 In g-^(5(x2 — X2) 



Xn 
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5{xi — x\] 



'\n{l- z) 



5{xi — x\) 



X2 



'ln(l — a;2/x2) 







X2 -XI 



(A.5) 



+a;° ("li^lx") - ^ + 1 111^(1 - xl)\ 6{X2 - 4) 



and similar for (1 ^^ 2). These identities are however inapphcable to the only remaining 
term 



l-z) 



(1 - z){xi + X2) 



{ztU (1-4) + . 

as ln(xi — Xi) and ln(x2 — x^) appearing in ( [A.4|) are not defined for z = 1. To overcome 
this problem, we rewrite 



ln(xi — x^] 



^0 
X2 
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1 -^1 X]^X2 / J 
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and similar for (1 ^^ 2). After subsequent algebraic manipulations and integration over 
ai and 0^2, we obtain: 
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(A.6) 



Analogous relations can be obtained for the ^-distributions, where we can first make 
the replacement 



6{l-z)6{y 



j^ a^ j^ ^ + 2:1 - zzi 



X2 



1 — Zi+ ZZi 



x5x25(Xi — X^)6{X2 — X2). 



Evaluation of the ^i-integral in the remaining terms yields 
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'^ ^2(1 — z)(XiX2 + X1X2) ^ ^ Xi(l — z)(XiX2 + X1X2) 

Only the identities (|A.5|) are unchanged for the ^-distributions, the remaining identities 
read: 
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